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ABSTRACT 

<N 

Aims. We aim to study the physical conditions, particularly the excitation state, of the intermediate-temperature gas in protoplanetary nebulae and 
young planetary nebulae (PPNe, PNe). The information that the observations of the different components deliver is of particular importance for 
understanding the evolution of these objects. 

Methods. We performed Herschel/HIFI observations of intermediate-excitation molecular lines in the far-infrared/submillimeter range in a sample 

00 ' of ten nebulae. The high spectral resolution provided by HIFI allows the accurate measurement of the line profiles. The dynamics and evolution 
£Nj , of these nebulae are known to result from the presence of several gas components, notably fast bipolar outflows and slow shells (that often are the 

fossil AGB shells), and the interaction between them. Because of the diverse kinematic properties of the different components, their emissions can 
• 1 be identified in the line profiles. The observation of these high-energy transitions allows an accurate study of the excitation conditions, particularly 
in the warm gas, which cannot be properly studied from the low-energy lines. 
, Results. We have detected FIR/sub-mm lines of several molecules, in particular of 12 CO, 13 CO, and H 2 0. Emission from other species, like 
NH3, OH, H2 18 0, HCN, SiO, etc, has been also detected. Wide profiles showing sometimes spectacular line wings have been found. We have 
mainly studied the excitation properties of the high-velocity emission, which is known to come from fast bipolar outflows. From comparison with 
general theoretical predictions, we find that CRL618 shows a particularly warm fast wind, with characteristic kinetic temperature 7k ; 200 K. In 
^ , contrast, the fast winds in OH 23 1 .8+4.2 and NGC 6302 are cold, T k ~ 30 K. Other nebulae, like CRL 2688, show intermediate temperatures, with 
5^ ■ characteristic values around 100 K. We also discuss how the complex structure of the nebulae can affect our estimates, considering two-component 
■ models. We argue that the differences in temperature in the different nebulae can be due to cooling after the gas acceleration (that is probably due 
' to shocks); for instance, CRL 618 is a case of very recent acceleration, less than ~ 100 yr ago, while the fast gas in OH 23 1.8+4.2 was accelerated 

1 1 , ~ 1000 yr ago. We also find indications that the densest gas tends to be cooler, which may be explained by the expected increase of the radiative 

. cooling efficiency with the density. 

^ . Key words, stars: AGB and post-AGB - stars: circumstellar matter, mass-loss - planetary nebulae 
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1 . Introduction during the AGB phase (e.g. lBalick & Frankl2002l) . The presently 

observed bipolar outflows would then correspond to a part of the 



Q\ ■ Asymptotic giant branch (AGB) stars, together with the thick cir- re i a tively dense shells ejected during the last AGB phase, mostly 
cumstellar envelopes (CSEs) around them, evolve to form plan- their polar re g ionSi accelerated by the shocks that propagate dur- 
i etary nebulae (PNe), through the phase of protoplanetary neb- m the PPN hase M ects of this ss ^ however 
t-H , ulae (PPNe). AGB stars are mass-losing red giants surrounded no j y e j we jj understood 
!LJ ■ by spherical, slowly expanding CSEs, while PNe, surrounding 
. ^ blue or white dwarfs, show in most cases strong axial symmetry Tne massive bipolar outflows in PPNe, as well as the unal- 
^ . or ring-like geometries and very fast bipolar jets. The transition tere d remnants of the AGB shells, usually sho w strong emission 
H ; time between these very different stages is extremely short, only m mo lecular lines (e.g. iBuiarrabal et alEoOlh . which is the best 
, CS , about 1000 yr, during which the intermediate PPNe develop a too i to stu d y these nebulae. Several PPNe have been extensively 
wide variety of structures. observed in mm-wave lines, including interferometric maps with 
PPNe and young PNe already present fast bipolar outflows, resolutions Thanks to those observations, the structure, dy- 
along with slower components with velocities similar to those of namics, and physical conditions in a number of PPNe are quite 
AGB envelopes. These bipolar flows typically reach velocities of well known. The profiles of the molecular lines have been found 
50-100 km s~' , and affect a sizable fr action of the nebular mass, to be composite, showing in most sources a central core plus 
~ 0.1-0.3 M e (Buiarrab al et al.ll200ll) . Such remarkable dynam- wide line wings. Maps indicate that the line core basically comes 
ics is thought to be the result of the interaction between the AGB from the remnant of the slow AGB wind, which is found to be 
and post-AGB winds: axial, very fast jets ejected in the very first roughly spherical or confined to equatorial regions. On the other 
post-AGB phases would collide with the denser material ejected hand, the line wings come from the fast bipolar outflows, which 
are more or less elongated in the direction of the nebular axis 



* Herschel is an ESA space observatory with science instruments pro- and show in most sources a strong velocity gradient, the ve- 
vided by European-led Principal Investigator consortia and with impor- locity incre asing proportionally to the distance to the s t ar. Se e 
tant participation from NASA. examples in Sanc hez Contreras et al.l (120041) . iFong et a"D (j2006), 
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lAlcolea et all (120071) . ICastro-Carrizo et al.l (I2010I) . and further 
references given below. 

However, observations of the low-/ transitions are not very 
useful for studying the warm gas components, in particular those 
heated by the shocks. The often observed J =2-1 and 7=1-0 
transitions of CO only require temperatures 7\ ~ 15 K to be 
excited. Indeed their maximum emissivity occurs for excitation 
temperatures of 10-20 K, and the line intensities and line inten- 
sity ratios depend only slightly on the excitation state in rela- 
tively warm gas. Of course, observations in the visible or near 
infrared ranges tend to select hot regions, with typical tempera- 
tures over 1000 K. The proper study of warm gas, 100 K < 7\ 
5 1000 K, therefore requires observations at intermediate wave- 
lengths, in the far infrared (FIR) and sub-mm ranges. 

Because of the crucial role of shocks in the formation of PNe, 
these warm regions are particularly important for understand- 
ing the nebular structure and evolution. In so me well-studied 
cases e.g. M 1-92, OH 23 1 .8+4.2. and M 2-56 dBuiarrabal et al l 
119981 lAlcolea et al.ll200ll ICastro-Carrizo et alJl2002h . the high 
velocity, massive outflows seem to be very cold, with tempera- 
tures < 20-30 K, which implies very fast cooling in the shock- 
accelerated gas. No warm component, presumably heated by re- 
cent shocks, has been identified in these sources until now. 

The only PPN in which detailed studies of warm acceler- 
ated gas have been peformed is CRL618. Interfero metric imag- 
ing of the I2 CO 7=2-1 line shows that dense gas in axial 
structures presents higher tem peratures, typically ~ 50-100 K 
(San chez Contreras et al.1120041) . But, precisely because of their 
relatively high excitation, the temperature estimate in these com- 
ponents from CO 7=2-1 is uncertain. Maps of 12 CO 7=6-5 
were also obtained by Nakashima et al. (2007), but with much 
less information on the nebular structure. A first attempt to study 
such w arm regions in the FIR was performed bv lJusttanont et al.l 
(2000), from ISO measurements of CO lines as high as 7=37- 
36 in CRL 618, CRL 2688, and NGC 7027. These authors indeed 
found temperatures of several hundred degrees in these sources. 
The lack of spectral resolution in the ISO data prevented any 
discussion on the components responsible for the line emission, 
which made difficult the interpretation of these results in nebu- 
lae showing very complex structure and dynamics. (This work 
and other previous ones are discussed in more detail in Sect. 
3.2.) Rece nt Herschel/HIFI obser vations in the FIR/sub-mm of 
CRL 618 dBuiarrabal etal.ll2010l) . with sufficiently high spec- 
tral resolution, have confirmed the presence of high temperature 
in various kinematic compon ents. The temperature derived by 
San chez Contreras et al ] (120041) was found to be too low, the fast 
gas showing temperatures typically ~ 200 K in the bipolar out- 
flow, and significantly higher in its base. Similar high tempera- 
tures were found in the innermost regions of the slow compo- 
nent. 

The Herschel Space Telescope is well-suited for studying 
warm gas around evolved stars in the FIR and sub-mm. The high 
spectral resolution that can be achieved with the heterodyne in- 
strument HIFI (better than 1 km s _1 ) is particularly important in 
our case: it is useful to identify the different nebular components 
in the profiles and allows the analysis of the kinematics, a funda- 
mental parameter in the study of this warm, shocked gas. Here 
we present Herschel/HIFI observations in a sample of ten PPNe 
and young PNe of several intermediate-excitation lines of abun- 
dant molecules: I2 CO and I3 CO, H 2 0, OH, NH 3 , HCN, SiO, etc. 
These observations have been obtained as part of the Herschel 
guaranteed-time key program HIFISTARS, which is devoted to 
the study of intermediate-excitation molecular lines in nebulae 
around evolved stars. 




1151 905 c,H/. i:;h 




o91.1?:.l GHz USB 




Fig. 1. HIFI observations including detected 12 CO, 13 CO, and 
C ls O lines in CRL 618 (7 m b vs. LSR velocities). The shadowed 
regions indicate the velocity ranges in which we can expect 
emission from the considered line, often mixed with emission 
from other lines. Note for instance the emission of HCN super- 
posed to 12 CO 10-9 and the overlapping lines 12 CO 16-15 and 
OH 3/2-1/2 (for more details, see the observed full-band spectra 
in the Appendix). 



2. Observations 

We used the HIFI heter odyne spectrometer, on board the 
Hersc hel Space Telescope (Pilbra tt et alJl2O10L Ide Graauw et al.l 
2010), to observe a number of molecular lines in a total of ten 
protoplanetary nebulae (PPNe). The observed objects are given 
in Table 1, where we quote the observed coordinates and addi- 
tional properties of the sources. 

Herschel/HIFI is an instrument very well suited to observe 
molecular lines in the FIR and sub-mm domains, providing a 
very high spectral resolution that allows the line profiles to be 
resolved. The data were taken using the two orthogonal HIFI 
receivers, H and V, available at each band, which were systemat- 
ically averaged after checking that there was no significant dif- 
ference in the lines detected in both receivers. The spectrometer 
works in double side-band (DSB) mode, which effectively dou- 
bles the instantaneous IF coverage. Care was taken when choos- 
ing the local oscillator frequency, maximizing the number of in- 
teresting lines that could be observed simultaneously. 

A total of nine frequency settings were used (not all of them 
in all sources), between 557 and 1843 GHz. The settings were 



V. Bujarrabal et al.: Herschel/HIFI observations of molecular emission in protoplanetary nebulae and young planetary nebulae 
Table 1. Observed sources. 



3 



source 


a (J2000) 


S (J2000) 


central star 


previous CO data 


comments; distance 


C-rich nebulae 


CRL618 


04 42 53.7 


+36 06 53.2 


B0 


2-1 and 1-0 maps 1 , high- J data 2 




CRL2688 


21 02 18.7 


+36 41 37.8 


F5Iae 


2-1 maps 3 




NGC7027 


21 07 01.6 


+42 14 10.2 


O 


2-1 and 1-0 maps 4 


high-excitation nebula 


Red Rectangle 


06 19 58.2 


-10 38 14.7 


Al 


2-1 and 1-0 maps 5 


equatorial disk in rotation 


IRAS 22272+5435 


22 29 10.4 


+54 51 06.3 


G5Ia 


1-0 maps 6 , 2-1 and 1-0 spectra 7 


isotrpic shell, no bipolar flow 


O-rich nebulae 


OH 23 1.8+4.2 


07 42 16.8 


-14 42 52.1 


M9III + AV 


2-1 and 1-0 maps 8 , high-./ data 9 


Calabash Nebula 


Boomerang Nebula 


12 44 45.4 


-54 31 11.4 


G0III 


2-1 and 1-0 spectra 10 




Frosty Leo 


9 39 54.0 


+ 11 58 54.0 


K7III 


2-1 and 1-0 maps 11 




IRAS 17436+5003 


17 44 55.4 


+50 02 39.5 


F2-5Ib 


2-1 and 1-0 spectra 7 


HD 161796 


NGC 6302 


17 13 44.2 


-37 06 15.9 


O 


2-1 maps 12 


high-excitation nebula 



Refer ence s: ': ISanchez Contreras et all d2004|): 2 : IBuiarrabal et alj i201Ch. INakashima et alj d2007h: 3 : 
J2010T): 5 : IBuiarrabal et al.1 d2005h: 6 : Fong et al] d2006h: 7 : iBuiarrabal et alj d2001h : 8 : lAlcolea eTal. ._ 
Buiarrabal & Bachillerl dl99lh . lSahai & Nvmanl dl997h : ' 1 : ICastro-Carrizo et aiTi2005h : 12 : lDinh-V-Trung et al.1 12008) 



Cox et al] d200Ch: 4 : INakashima et all 
d2001h; 9 : lAlcolea etal] d2011h : ' 



CHL 618 H a O, OH & NH 3 
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JSB fl 






OH 2 n V3 .1=3/2 -l/2 + 


1834.75 GHzfll 
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1163.912 GHz 
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h 1113.343 GHz USH 5] 




o-H 3 J^,^=li,o-lo,j 
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Fig. 2. HIFI observations of CRL618 including detected lines of 
H2O, OH, and NH3 (T m \, vs. LSR velocities). The shadowed re- 
gions indicate the velocity ranges in which we can expect emis- 
sion from the considered line, often mixed with emission from 
other lines (see identifications in Fig. A.l). 

chosen to observe the 12 CO and 13 CO 7=6-5, 7=10-9, and 
7=16-15 lines and the H 2 li,o-lo,i. li,i-Oo,o, 3 2 ,i-3i, 2 , 3i, 2 - 
2 2 ,i, 42,2— 4i,3, and 73 4-72,5 lines. Due to the wide simultaneous 
frequency bands observable with HIFI and the double side band 
system used, each frequency setting often include other lines. 
Interesting molecular transitions also placed within the observed 
bands are the NH 3 7,tf=l,0-0,0andOH 2 ni /2 3/2-1/2 rotational 
lines, as well as lines of C 18 0, H 2 ls O, vibrationally excited H 2 0, 
SiO, HCN, CN, etc. 



CRL 2688 12 C0, 13 C0 & C 18 



' ld CO 7=16-15 




1760-486 GHz USB 


T-^rwn" ^ n^y w ^ »[ 


1E CO 7=16-15 




1841.346 GHz USB 




L3 C0 7=10-9 




1101.350 GHz LSB 




I I3 C0 7=10-9 




1151.985 GHz LSB ; 




7 C 18 7=6-8 




658.506 GHz USri : 


13 go j=6-5 r^y\ 


661.067 GHz USB " 




1 la C0 7=6-5 




691.473 GHz USB J 



-150 -100 -00 50 

V,. SR (km s ) 



Fig. 3. HIFI observations including detected 12 CO, 13 CO, and 
C ls O lines in CRL 2688 (T mh vs. LSR velocities). 



The observed full-band spectra are reproduced in the 
Appendix. In those figures we indicate the upper and lower side 
band (USB and LSB) frequency scales and the frequencies of in- 
teresting lines placed within the observed bands. Continuum has 
been subtracted in all the spectra presented in this paper. 

The observations were obtained in the dual-beam-switching 
(DBS) mode. In this mode, the HIFI internal steering mirror 
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NGC 7027 12 CO, 13 C0 & C 1B 

Vb (K) 




-50 _ l 50 100 

Vlse (km s ) 



Fig. 4. HIFI observations including detected 12 CO, 13 CO, and 
C ls O lines in NGC 7027 (T mb vs. LSR velocities). 
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- ' OH E n 1/2 J=3/2 


j| n rLnn 1834.75 GHz LSB 

ft \t*%. 






I 














; p-H 2 J" Ka ,K c =l u -0 0>0 rH 11 13.343 GHz USB 













-50 50 100 

(km s" 1 ) 



Fig. 5. HIFI observations of NGC 7027 including detected lines 
of H2O and OH (T m \, vs. LSR velocities). 



chops between the source position and an emission-free position 
3' away. The telescope then alternatively locates the source in 
either of the chopped beams, providing a double-difference cal- 
ibration scheme, which allows a more efficient cancellation of 
the residual standing waves in the spectra. This procedure works 
very well except for the highest frequencies (HIFI band 7, above 
1700 GHz), where noticeable ripples were found in some spec- 
tra, specially in the V-receiver. Other spurious effects in the base- 
lines were found in a few cases, particularly around 1200 GHz 
due to very different behaviors of the sub-band baselines, which 
were sometimes difficult to remove. 



Red Rectangle CO & CO 

T mt (mK) 




-40 -20 20 40 

V LSR (km s ) 



Fig. 6. HIFI observations including detected CO and CO 
lines in the Red Rectangle (r m b vs. LSR velocities). 



IRAS 22272 + 5435 12 C0 & ,8 C0 

T„b (mK) 




-100 -50 50 

v mi< (km s ) 



Fig. 7. HIFI observations including detected 12 CO and 13 CO 
lines in IRAS 22272+5435 (r mb vs. LSR velocities). 



The HIFI data shown here were taken using the Wide-Band 
Spectrometer (WBS), an acousto-optical spectrometer that pro- 
vides simultaneous coverage of the full instantaneous IF band 
in the two available orthogonal receivers, with a spectral reso- 
lution of about 1.1 MHz. In some cases, the spectra have been 
smoothed to improve the S/N ratio. 
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-100 100 20C 

Va, (km s ) 



Fig. 8. HIFI observations including detected 12 CO, 13 CO, and 
C ls O lines in OH 23 1.8+4.2 (r mb vs. LSR velocities). Note that 
other lines appear in some spectra, see Fig. A. 6. 



We processed the data using the standard HIFI pipeline in the 
HIPE software, with a modified version of the level 2 algorithm 
that yields unaveraged spectra with all spectrometer sub-bands 
stitched together. Later on, the spectra were exported to CLASS 
using the hiClass tool within HIPE, for further inspection, "flag- 
ging" data with outstanding ripple residuals, final averaging and 
calibration, and baseline removal. Even if spurious effects in the 
baselines were in a few cases impossible to remove, we have 
only considered baselines of moderate degree, in general just of 
degree 1, at most of degree 3. 

The data were originally calibrated in antenna temperature 
units and later converted into main-beam temperatures (T m \,). 
A summary of the telescope characteristics and observational 
uncertainties, including the spatial resolution and the conver- 
sion factors to main-beam units, can be found in the Herschel 
Observers' Manual and the HIFI Instrument and Calibration 
Web Page (see http://herschel.esac.esa.int). We recall that the 
FWHM beam sizes range between about 33 and 12 arcsec (in- 
versely proportional to the observed frequency) and T^/T^ ra- 
tios range between 1.3 and 1.5. 

In Figs. 1 to 14, we show zoom-ups of selected lines de- 
tected in our sources, in LSR velocity and T m \, units. In CRL 618, 
some of these l i nes w ere already presented and discussed by 
Bujarraba l et alJ ([2010); here we have reanalyzed all the data 
and included new observations. In Figs. 1 and 2, corresponding 
to CRL 6 18, we have indicated in grey the velocity range ex- 
pected for each line profile, because of the very complex spectra 
and wide profiles characteristic of this nebula. Data of the most 
interesting sources will be described and analyzed in detail in 
forthcoming papers. 




-100 100 20C 

V m (km s -1 ) 



Fig. 9. HIFI observations of OH 23 1.8+4.2 including detected 
lines of H2O (T m b vs. LSR velocities). Note that other lines ap- 
pear in some spectra, see Fig. A. 6. 



Boomerang Nebula CO & H 2 

Tab M 




-100 100 

V 1SH (km s >) 



Fig. 10. HIFI observations including detected CO and H2O 
lines in the Boomerang Nebula (T m \, vs. LSR velocities). 



3. Results 

In several objects, namely in CRL 6 18, CRL 2688, 
OH 23 1.8+4.2, NGC6302, Frosty Leo, and IRAS 17436+5003 
(HD 161796), we can separate the emission of the line central 
core (which is the dominant peak in almost all cases) from that 
of the line wings, which are expected to come, respectively, 
from the fossil AGB envelope and from the fast outflows result- 
ing from the post-AGB axial acceleration (S ect. 1). We have 
alread y studied in detail the case of CRL 618 (Bujarraba l et al.l 
120101) : our data of CRL 2688 and OH 23 1 .8+4.2, which are also 
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Frosty Leo 12 C0, 13 C0 & H 2 

T mb (mK) 




-200 -100 , 100 

V„ (km a" 1 ) 



Fig. 11. HIFI observations including detected 12 CO, 13 CO, and 
H2O lines in Frosty Leo (T m b vs. LSR velocities). 



IRAS 17436 + 5003 12 C0 & 13 C0 

T mb (mK) 




-150 -100 -50 50 100 

V^p (km a l ) 



Fig. 12. HIFI observations including detected I2 CO and 
13 CO lines, r mb vs. LSR velocities, in IRAS 17436+5003 
(HD 161796). 



well known nebulae, will also be analyzed more in depth in 
forthcoming papers. 

In others (the Red Rectangle, the Boomerang Nebula, 
and IRAS 22272+5435), no line wings are detected. The 
Red Rectangle shows a peculiar dynamics dominated by 
rotati on, involving only moderate speeds dBuiarrabal et al.l 
120051) : see preliminary analysis of our data in Sect. 3.1. In 
IRAS 22272+5435, no fast flows have been detected to date 
and molecules are apparently abundant only in a relativel y 
slow-moving shell (Buiarrabal et~ai1 l200ll iFong et al.l 12006). 
The Boomerang Nebula, a less well studied object, shows a 
low-velocity warm component plus, probably, very cold high- 
velocity gas, which cannot be detected in high- 7 CO emission 
dBuiarrabal & Bachilleri[l99"ll ISahai &Nvmanll 19971) . 

Finally, NGC7027 is a complex nebula mainly composed 
of a s lightly elongated (empt y) shell that shows many peculiar- 
ities (Nakashima et al. 2010). Our line profiles are composite 
and intense, as expected in view of mm-wavelength data. In this 
source, we did not try any simplified discussion in terms of line- 
core and -wing emission because of its complex structure. The 
interpretation is also hindered by the relatively large angular ex- 
tent of the nebula. Our single-point observations of CO J -2-1 
(from the 30m telescope) and 7=16-15 (from Herschel), those 



NGC 6302 12 CO & 13 C0 

T„ (K) 




-100 -50 

V 1SH (km sT 1 ) 



Fig. 13. HIFI observations including detected 12 CO and I3 CO 
lines in NGC 6302 {T m b vs. LSR velocities). 
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Fig. 14. HIFI observations including detected lines of H2O, OH, 
and NH 3 in NGC 6302 (r mb vs. LSR velocities). 

showing the highest angular resolution, are more sensitive to the 
central part of the nebula image. Therefore, these data tend to de- 
tect only the extreme LSR velocities (~ 10 and 40 km s~'), since, 
due to projection effects, the emission in these velocities comes 
from the central regions of the expanding shell image. The op- 
tically thin 13 CO lines and C ls O 7=6-5 show high intensity in 
the central velocities; they should come from slow gas with very 
high column density. 12 CO 7=6-5 shows a very strong self ab- 
sorption in the blue part of the spectrum, expected in optically 
thick lines from expanding gas. This self absorption is much 
stronger in this line than in the 7=2-1 and 7=1-0 ones, indicat- 
ing the presence of particularly high opacity and relatively warm 
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gas in the shell. Finally we note that H2O lines clearly come 
from the fastest gas, though in this case the spatial resolution is 
not very high. These results suggest that H2O formation in this 
C-rich environment is a result of the shocks that accelerated the 
gas. Curiously, the OH emission (observed with higher spatial 
resolution) is not intense in high-velocity features. Obviously, 
NGC 7027 is one of the sources for which the high-excitation 
molecular emi ssion requires a careful mode ling; results will be 
soon reported (Sant ander-Garcfa et alJl20l"Tl in preparation). 

In Table 2, we show characteristic intensities of the CO lines 
averaged over representative regions of the profiles, the center 
and the blue and red wings, in the sources in which they can be 
well defined. Intensities of H2O lines in the same velocity ranges 
are shown in Table 3. The adopted velocity ranges are given in 
those tables. In order to follow a procedure as objective as pos- 
sible, avoiding any bias due to a priori ideas about the nature 
of the sources, intensities are averaged in all cases over velocity 
ranges 5 km s wide and defined following a standard proce- 
dure. This width seems a good compromise to yield an intensity 
average over a relatively large number of spectral channels, but 
being narrow enough to allow the separation of each spectral fea- 
ture, as we will see below, and to avoid integrating spectral re- 
gions with too weak emission. In the definition of these velocity 
ranges, we have taken into account previous observations of low - 
7 CO lines to identify the different components in the line pro- 
files, which are not always well seen in our profiles. In particular, 
values of the systemic and expansion velocities of the slow com- 
ponents ar e mainly based on d ata and estimates of these param- 
eters by Bujarrabal et afl (l200lh . We have also confirmed from 
published maps the identification of such features in terms of 
kinematical components and checked possible contaminations 
between their emissions. For the line wings, the representative 
velocity slots were slightly more separated from the center than 
the slow-component expansion velocity, to avoid contamination 
from that component, in particular the relative minima due to ab- 
sorption by the outer slow gas, which can be clearly seen in some 
profiles. In a few cases, the line core may also be contaminated 
by emissions of fast outflows that, due to projection effects, may 
show LSR velocities close to the systemic one; we have checked 
for each individual nebula, taking into account the existing in- 
formation (see below), that this contamination is negligible for 
the used velocity slots. 

We can see detaile d information on the structure and kine- 
matics of CRL618 in Sanch ez~Contreras et al.l (120041): maps o f 
IRAS 17436+5003 are shown by ICastro-Carrizo et alj ((2004). 
Both sources show an extended slow component, inside which 
fast bipolar flows are placed, whose emissions can be well iden- 
tified in single-dish profiles (Bujarrab al et al.l [2.001b . The gen- 
eral structure of CRL2688 is similar, but the fast gas struc- 
ture is com plex and several c ollimated jets exist, see lFong et alJ 
(200fi|) and dCox et a l. 2000, note that the extended component 
is significantly filtered out in these interferometric data). The 
relatively narrow velocity slot in the center used here mini- 
mizes contamination by emission from jets running almost in 
the plane of the sky, which can be present in this source; on 
the other hand, the contribution of the fast component to the 
single-dish profiles is well delimited by the absorption feature 
seen in single-dish observations. The fast bipolar flows in Frosty 
Leo are more extended than the slow molecular gas (though 
smaller than the optical nebula), emissions fro m both compo- 
nents being clearly separated in L SR velocity (Bujarrabal et alJ 
1200 1L Cast ro-Carrizo et alJ 120051) . The resulting LSR velocity 
ranges taken for CRL618 are [-24.5,-19.5] kms 1 (line cen- 
ter), [-49,-44] km s -1 (blue wing) and [0,+5] km s _1 (red 



wing); for CRL2688, we took respectively [-37.5,-32.5], 
[-62,-57], and [-13,-8] kms -1 ; for Frosty Leo we took re- 
spectively [-13.5,-8.5], [-27,-22], and [0,-1-5] kms" 1 ; and 
for IRAS 17436+5003 we assumed respectively the ranges 
[-37.5,-32.5], [-54,-49], and [-21,-16] km s _1 . 

The case of OH 23 1.8+4. 2 is different, because it shows a 
very elongated CO-rich component w ith a very well-defined ve- 
locity gradient dAlcolea et al.l l200"lh . in which the fossil AGB 
envelope is not clearly identified. In this nebula, we selected 
the LSR ranges [+30.5, +35.5] (line core), [0,+5] (blue wing), 
[+61, +66] (red wing). The two high-velocity ranges correspond 
to regions clearly different from the cent ral condensation (in 
clumps 12 and 14 in the nomenclature of lAlcolea et al.l [2001 . 
which have been very probably accelerated along the nebular 
axis), but not very far from the center, to be sure that their 
emissions are intense and well within the telescope beamsize. 
The central velocity range we selected represents the central 
clump, which, in view of its structure, kinematics, and chemistry, 
does not seem to be affected by pos t- AGB shocks (lAlcolea et al.l 
200lL lSanchez Contreras et al 1 fl997h . 

In NGC 6302, most of the molecular gas occupies an equa- 
torial torus or disk, in moderate expansion (~ 15 kms~') and 
perpendicular to the spectacu l ar opt ical lobes, see maps in 
mm-w ave CO by iPeretto et al.l J2007) and Dinh-V-Trung et all 
(2008). These authors argue that the equatorial structure is prob- 
ably the remnant of the (late) AGB mass loss, not strongly al- 
tered during the post- AGB phase. The two-horn dominant fea- 
ture of our profiles appears as a result of the shape and rela- 
tively constant velocity of this component. There are also high- 
velocity clumps extending partially along the lobes; but they are 
very patchy and only those at relatively negative velocity (as- 
sociated to the southeast part of the nebula) are intense enough 
to be analyzed here. Efficient photodissociation due to the UV 
radiation from the very hot central star is probably the reason 
for the distribution of molecule-rich gas found in NGC 6302. 
Therefore, in this source we just considered its blue wing, range 
[-53,-48] km s _1 , and the line core, range [-25,-20] km s _1 . 
We will assume that these velocity ranges represent respectively 
the high-velocity molecular gas and the unaccelerated remnant 
of the AGB shell, but keeping in mind that the nature of both 
components in this source may be different than in other nebulae 
studied here and the existence of large amounts of nebular gas 
not detected in molecular emission. 

We have computed wing/core CO intensity ratios and ra- 
tios between the observed lines and previous data of the 
CO 7-2-1 transitions (for previous observations, see Table 1, 
Bujarraba l et alJ200ll and references therein). For the wings, we 
have taken the average of both line wings, except for NGC 6302, 
in which only the blue wing is considered. See results in Table 2 
and in Figs. 15, 16, 17. It is obvious from these ratios that there 
is a wide variety of excitation conditions in the high-velocity 
(presumably shocked) gas. In Fig. 15 we can see the wing/core 
intensity ratios for three well known sources. These variations 
are not due to significant differences in the central core intensi- 
ties, for which the ratio with respect to the 2-1 transition varies 
very slightly from line to line and from source to source. The 
12 CO and 13 CO line-wing ratios with respect to the 7=2-1 line 
are shown in Fig. 16, these ratios are also a measurement of 
the excitation state of the high-velocity gas. The fast increase 
of the wing/core ratio with the energy of the levels in CRL618 
shows that the temperature in its fast outflow must be high; the 
case of OH 23 1.8+4. 2 is completely different, with a clear de- 
crease of this ratio, showing a cool fast flow; CRL2688 seems 
to be an intermediate case. For CRL618 and OH231.8+4.2 we 
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Table 2. Characteristic intensities of the central core (or peak) and wings of CO lines in the sources in which they can easily be 
identified. 
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Notes: See discussion on the respective velocity ranges in the text. Ratios wing/core and with respect to the 7=2- 1 lines (from 30-m telescope 
data) are also indicated; in these cases, the average of both wings is taken. 



have an accurate idea of the temperature in the outflows from de- 
tailed model fitting of Herschel CO data, which indeed confirms 
this preliminary interpretation. Bujarrabal et al. (2010) showed 
that the temperature of the fast flow of CRL 618 is high, 7\ 
Z 200 K. However, we know that the fast flow t emperature in 
OH 2 3 1.8+4.2 is always cooler than about 30 K dAlcolea et alj 
l201ll in preparation). 



3.1. Comparison with simplified calculations of CO line 
emission 



Our line profiles must be compared with theoretical predictions 
in order to derive the physical properties of the nebular gas from 
the empirical data. Here we just present a simplified descrip- 
tion of the analysis to be performed, mostly based on single- 
component, standard calculations. More detailed modeling of 
the nebulae will be performed in dedicated papers, where their 
complex structures will be taken into account. 
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Table 3. Characteristic intensities of the central core (or peak) and wings of H2O and OH lines in the cases in which they can be 
measured. 
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Fig. 15. Line-wing over line-core main-beam temperature ratio 
of the observed CO lines in the well studied sources CRL 6 18, 
CRL 2688, and OH 23 1.8+4.2 (7=2-1 data are taken from pre- 
vious works); see Table 2. Note the different behavior of the rel- 
ative line- wing intensity between CRL 618, which shows partic- 
ularly hot high-velocity outflows, and OH 23 1.8+4. 2, in which 
the flows remain quite cold. CRL 2688 is an intermediate case. 
Squares (continuous line) represent 12 CO data and triangles 
(dashed lines) show 13 CO data. Each rotational transition is rep- 
resented by the energy (in K) of its upper level. 



A comparison of the wing intensity ratio between our lines 
and 7=2-1 (see Fig. 16) with theoretical predictions can be used 
to perform a rough estimate of the physical conditions in the fast 
outflows in our sources, particularly the characteristic tempera- 
ture. Results from plain rotational diagrams, which require op- 
tically thin and thermalized lines, are not accurate in our cases, 
since lines a high as 7=16-15 are only thermalized for densities 
J 10 7 cirT 3 and optically thin emission is not expected in 12 CO 



200 400 600 line exc. (K) 

Fig. 16. Wing intensity ratio between the observed lines and the 
CO 7=2-1 transition (from previous works) in the well stud- 
ied sources CRL618, CRL 2688, and OH231. 8+4.2; see Table 
2. A comparison between these ratios and estimates from our 
simple, general models reveals values of the characteristic tem- 
perature compatible with estimates from detailed models (see 
text). Squares (continuous line) represent 12 CO data and trian- 
gles (dashed lines) show 13 CO data. Each rotational transition is 
represented by the energy (in K) of its upper level. 



lines (even 13 CO lines seem often to present moderate values of 
the optical depth). 

To perform a very simplified estimate of the excitation condi- 
tions in our nebulae, we have carried out theoretical calculations 
of line ratios for a variety of values of the density and tempera- 
ture (n tot and 7\), using the well known LVG (large velocity gra- 
dient) formalism. We have assumed standard values for the rest 
of the parameters that enter the calculation: distance to the center 
R = 6 10 16 cm, expansion velocity V e xp = 30 km s _1 , and relative 
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Fig. 17. Line-wing over line-core main-beam temperature ra- 
tio of the observed CO lines in sources that are less well 
studied in molecular emission: Frosty Leo, NGC 6302, and 
IRAS 17436+5003. Also in these cases, we can identify hot- 
ter and cooler fast outflows. Squares (continuous line) represent 
12 CO data and triangles (dashed lines) show 13 CO data. Each ro- 
tational transition is represented by the energy (in K) of its upper 
level. 



abundances X( 12 CO) = 2 10" 4 , X( 13 CO )= 10~ 5 . It is obvious 
that those conditions cannot be applied to all nebulae, since we 
assume a single value for the product R X/V exp (n tot R X/V exp rep- 
resents the column density in the opacity under the LVG approx- 
imation). However, this approach is enough for our purpose of 
understanding the general properties of the line intensities from 
the theoretical point of view, particularly to give examples of the 
behavior of the line intensities and ratios depending on the tem- 
perature and density. Opacity effects are taken into account in 
the calculations, under the usual assumptions of this formalism. 

The LVG approximation itself is fully justified in our case, 
since it is well known that PPNe often show strong velocity gra- 
dients, with radial expansion velocities increasing more or less 
proportionally to the distance to the star. For this reason, we do 
not expect strong radiative interactions between points separated 
by a long distance and showing very different physical condi- 
tions. This approximation is appropriate for our purposes also 
because it does not require that details on the large-scale nebula 
structure are included in the code, provided that no radiative in- 
teraction between distant points takes place, allowing a general 
description of the line excitation for a number of different ob- 
jects. Moreover, in our bipolar flows, often elongated and show- 
ing large velocity gradients, the emission at each LSR velocity 
must come from a relatively compact region. In this way, the 
code gives characteristic values of the brightness for each line, 
from which we can readily calculate line intensity ratios. Line 
intensities at a given velocity can be also estimated if we assume 
an angular size for the emitting region. 

Results of our calculations for 13 CO and 12 CO are shown in 
Figs. 18, 19, 20, and 21. In our analysis, we will focus on the line 
intensity ratios, since ratios depend only slightly on the linear 
and angular extent of the source and are more easily described 
by general models. Figs. 18 and 19 show predicted line ratios 
for several densities and temperatures, within the ranges needed 
to explain our data. In this case we have assumed that the lines 
come from a region significantly smaller than the telescope spa- 
tial resolution (which is justified in our case), so that the com- 
parison of the predicted line brightness ratios must be corrected 
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Fig. 18. Predictions for the wing intensity ratios between the ob- 
served 12 CO lines and the 12 CO 7=2-1 transition obtained from 
our non-LTE calculations for our three lines, 7=6-5, 7=10-9, 
and 7=16-15, respectively dashed (black in electronic version), 
dotted (red), and dash-dot (green) lines. Intensity ratios are 
shown as functions of the kinetic temperature in the three panels 
for, respectively, densities equal to 10 7 , 10 6 



and 10 5 cm 3 . 



for the respective telescope beam widths (half-maximum beam 
widths of the Herschel observations are given in Sect. 2; we have 
adopted a value of 12" for the 30m-telescope observations of CO 
7=2-1). In Figs. 20 and 21, we show predicted intensities; nor- 
malized units are used to allow a better comparison between the 
intensities of the different lines and sets of physical conditions. 
The chosen units correspond to 7 m b(K) calculated assuming that 
the model nebula with density n tot = 10 6 crrT 3 occupies a region 
1 " wide, and taking the different telescope resolutions into ac- 
count (which are again much larger than the emitting region). To 
make easier the comparison between predictions for other den- 
sities, we have taken varying values of the emitting angular size, 
inversely proportional to the density, such that the total emitting 
mass is similar for the different model condensations. We can 
check that the intensities are of the order of those observed: for 
instance a model clump ~ 2" wide, with 7\ = 100 K and n tot = 
10 6 cm 3 would yield a typical brightness 7 mb ( 12 CO 7=10-9) ~ 
1 K, similar to some of our observed line-wing intensities. 

The dependence of the predicted high-7 line intensity on the 
temperature is clear, as expected, and can be seen both in line 
ratios (Figs. 18 and 19) and intensities (Figs. 20 and 21). The 
intensity of the relevant lines increases with temperature, but 
from a certain value of 7k it saturates and even decreases. Of 
course, the temperature limit of intensity saturation increases 
for higher values of 7, which imply higher level energies. We 
can also see that the line intensity ratios tend to increase with 
the volume density, particularly for the highest 7-values, but 
the dependence is less clear than with respect to the tempera- 
ture. The general properties of the excitation of the high-7 CO 
lines discussed here are therefore different than those of low- 
excitation transitions of high dipole moment species, like for in- 
stance HCO + , SiO, etc, often observed with ground-based tele- 
scopes. In these low-energy transitions with large radiative decay 
rates, the emission depends strongly on the density (the intensity 
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Fig. 20. 12 CO line intensities from our simplified model calcu- 
lations, for several values of the temperature and density, see 
intensity normalization procedure in Sect. 3.1. 



Fig. 19. Predictions for the wing intensity ratios between the ob- 
served 12 CO lines and the 13 CO 7=2-1 transition obtained from 
our non-LTE calculations for our three lines, 7=6-5, 7=10-9, 
and 7=16-15, respectively dashed (black in electronic version), 
dotted (red), and dash-dot (green) lines. Intensity ratios are 
shown as functions of the kinetic temperature in the three panels 
for, respectively, densities equal to 10 7 , 10 6 , and 10 5 cm -3 . 

per unit volume being roughly proportional to n 2 0t ), because col- 
lisional de-excitation is much less probable than radiative transi- 
tions. However, its dependence on the temperature is in fact very 
moderate, because the intensity saturation level is low. 

From comparison of our 13 CO observational results with 
model predictions in Fig. 19, we derive that the high-velocity 
CO emission in CRL618 essentially comes from a region at 7k 
~ 150 K and with a density n tot ~ 10 6 cm" 3 (as we will see below, 
our 12 CO data indicate that the temperature could be somewhat 
higher). In OH 23 1.8+4. 2, our data are compatible with a tem- 
perature ~ 30 K and a lower density. Both results are in good 
agr eement with detailed model estimates (see above discussion 
and lBuiarrabal et alJl20ldlAlcolea et al.ll201 lb . 

Our 12 CO observations of CRL618 are compatible with the 
predictions in Fig. 18 for a density of about 10 6 cm" 3 and a tem- 
perature of about 270 K, higher than the value deduced from 
13 CO data. This discrepancy may come from the fact that, in 
sources showing significant variations of the density within the 
region emitting at a given velocity, the optically thick 12 CO 
emission selects very dense regions less efficiently than the 13 CO 
lines, those denser regions being probably cooler (see discus- 
sion below on the line ratios predicted by two-component mod- 
els). [Bujarraba l et al.l (1201 Ob deduce from their detailed model 
fitting that the temperature of the high-velocity outflow must 
vary within this range, with typical val ues around 200 K, prob- 
ably h igher in the base of the outflow. San chez Contreras et al.l 
d2004b deduce a density similar to tha t given here (which was 
also assumed bv lBuiarrabal et al.ll20l6T) . For OH 23 1.8+4.2, we 
derive from the I2 CO lines 7\ ~ 30 K and n tot ~ 10 5 cm 3 , which 
are ag ain comparable to th e values obtained from detailed mod- 
els bv lAlcolea etaLl d2011h . 

In the case of CRL2688, 12 CO line ratios are compatible 
with a temperature of about 100 K and n tot ~ 10 6 cm" 3 . Such 
conditions are also roughly compatible with our 13 CO data. It is 



probable that also in this case the complex nebula harbors gas 
with variable physical conditions. 

As we see, even if the procedure depicted here is simple, 
the values we derive are quite self-consistent and consistent with 
other more accurate results from detailed modeling. 

We can also perform the same kind of analysis for 
the less well studied sources NGC6302, Frosty Leo, and 
IRAS 17436+5003. Obviously, fast winds in NGC6302 and 
Frosty Leo show a very low excitation. In Frosty Leo, 
ICastro-Carrizo et all d2005h predicted a low temperature, from 
low-7 observations, but our data indicate the presence of a 
slightly warmer component at high velocity, with a temperature 
~ 60 K (assuming n tot ~ 10 4 cm" 3 , following Castro-Carrizo et 
al.). In NGC 6302, the fast gas must present a low temperature ~ 
40 K, and the density must be of about a few 10 s cm" 3 ; both 
results are compatible with estimates by iDinh-V- Thing et all 
(12008b from interferometric 7=2-1 mapping. The fast flow in 
the less well studied IRAS 17436+5003 must be more excited, 
with 7k ~ 60 K and n tot ~ 10 6 cm" 3 . In this object, the properties 
of the molecular component and the evol utionary stage seem to 
be similar to those of CRL 2688 (see e.g. Bujarrabal et al.ll200H 
note that lines in IRAS 17436+5003 are weaker, mainly because 
it is more distant). 

As mentioned, these estimates are in any case crude, partic- 
ularly in view of the probably complex structure of the emitting 
regions. In order to investigate those effects in general terms, we 
must consider our calculations of intensity predictions in Figs. 
20 and 21. Let's take a two-component case. In general, we can 
see that when there is a massive component with « to t ~ 10 5 cm" 3 
and 7k 5 200 K the 7=2-1 emission is always dominant in it 
and, even if a denser/warmer component exists, the total inten- 
sity ratios 7=10-9/7=2-1 and 7=16-15/7=2-1 can hardly be 
larger than ~ 1/2. If there is a significant component with 7\ 
I 200 K and n M I 10 s cm" 3 , its 7=10-9 and 7=16-15 emis- 
sion would always be dominant over that of a diffuse/cool com- 
ponent. This would explain the general trend of relatively low 
7=10-9/7=2-1 and 7=16-15/7=2-1 ratios found in our data. 

We can also suppose that relatively diffuse, accelerated com- 
ponents are systematically warmer, if the cooling after shock 
passage is radiative (and, therefore, denser gas cools faster). This 
would yield regions with n tot ~ 10 5 cm" 3 and 7k J 200 K, whose 
emission is intense in 12 CO lines, but relatively weak in 13 CO 
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NORMALIZED LINE INTENSITY ( ,3 C0) vs. KINETIC TEMPERATURE 
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Fig. 21. CO line intensities from our simplified model calcu- 
lations, for several values of the temperature and density, see 
intensity normalization procedure in Sect. 3.1. 

lines (with respect to dense regions). After taking the emission 
of more cool/dense components into account, one would get rel- 
atively high 12 CO 7=10-9/7=2-1 and 7=16-15/7=2-1 line ra- 
tios, higher than the corresponding 13 CO ratios. Therefore, the 
presence of relatively weak line-wing emission in 13 CO high- 
7 lines, often found in our observational data, could be due to 
the presence of an important relatively diffuse component, de- 
tectable in 12 CO emission but practically not in I3 CO lines. 

The presence of selfabsorption by outer, less excited layers 
(not included in our simple discussion) could also help to ex- 
plain this trend. In this case, self absorption would affect mostly 
low-7 12 CO lines, since they present the highest opacities in low- 
excitation regions. This would yield relatively low I2 CO 7=2-1 
intensities and then relatively high 12 CO 7=10-9/7=2-1 and 
7=16-15/7=2-1 ratios. The existence of selfabsorption under 
these conditions is confirmed by the observed absorption fea- 
tures in mm-wave CO profiles. 

As mentioned, we have studied the case of CRL618 by 
means of detailed modeling of its line emission, based in our 
good knowledge of the stru cture of the nebulae from low-7 CO 
data (Buiarrab al et al.ll2010h . and we are analyzing in depth the 
high-7 lines of OH 23 1 .8+4.2, CRL 2688, and NGC 7027, which 
are also well known. 

We are also studying in detail the CO emission from the Red 
Rectangle, so here we just present a brief, simplified discussion. 
The Red Rectangle is the only PPNe or young PNe in which 
an equatorial disk in rotation around t he central star has bee n 
well detected (from mm-wave mapping, Buiarra bal et al.ll2~005l) . 
This disk is, as far as we know, the only molecule-rich com- 
ponent of the nebula. The structure of the line profiles in this 
source is accordingly peculiar, because the outer part of the pro- 
files come from the inner disk regions in fast rotation (rotation 
is exactly Keplerian in them), which are relatively hot and very 
small. Meanwhile, the central part of the profiles mainly comes 
from the outer disk. Let us take representative LSR velocities 
around the systemic one V± = ± 1.65 km s _1 (the LSR velocity 
of the center of gravity of this source is ~ km s ). These veloc- 
ities correspond to the rotation velocity at the point at which the 
dynamics change from pure Keplerian rotation to rotation plu s 
expansion (Rkep in the notation used by Buiarrabal et afl l2005). 
and there is a clear maximum of the CO lines in them. We can see 
in our spectra that 12 CO 7=6-5 and 7=10-9 show very similar 



values of the main-beam brightness at V±. This is true even for 
7=16-15, in spite of the poor-quality profile at this frequency. 
I3 CO 7=6-5 and 7=10-9 also show practically the same inten- 
sities at V ± . As we have discussed (Sect. 1), mm-wave CO data 
are adequate to measure densities and total masses, but often 
yield underestimates of the temperature because those lines tend 
to select the coolest compone nts. So, if we as s ume th e density 
estimated in these regions bv lBuiarrabal et alj d2005l) . 10 5 -10 6 
cm" 3 , we conclude from comparison with our predicted inten- 
sities, that the temperature in this representative region is > 150 
K, at least twice as high as the temperature deduced from mm- 
wave lines. This result has important consequences for the struc- 
ture and dynamics of the equatorial disk in the Red Rectangle. 
Under the framework of the standard theory of the passive ro- 
tating disks, the gas temperature directly controls the disk extent 
from the equator (but comparison with observations is difficult, 
because the angular resolution of the existing maps is not high 
enough to properly measure the disk width). The presence of 
these high temperatures is also the most probable explanation 
of the detection of evaporation of gas in this source (i.e. of ex- 
pansion from a certain radius R^ ep ); in fact, we can readily show 
that the thermal velocity dispersion becomes comparable to the 
rotation velocity at R^ ep for Ty_ ~ 200 K. 

3.2. Comparison with previous estimates 

As mentioned in Sect. 1, warm molecular gas in CRL 618 and 
CRL 2688 was stud i ed in a few previous papers. 

lJusttanont et al.l d2000t) deduce the presence of warm gas 
from ISO data of 12 CO emission, for transitions between 7=14— 
13 and 7=37-36. These authors clearly detect the presence of 
gas at temperatures higher than the typical values found here, of 
~ 380 K in CRL 2688 and ~ 700 K in CRL 618, using simple ro- 
tational diagrams. This result is not unexpected, because higher- 
7 transitions tend to come from hotter regions. However, the lack 
of information on the profile shapes (ISO spectra did not resolve 
the line profiles) prevents any attempt to discern from which neb- 
ular components such emissions arise. In the case of CRL 618, 
this hot component could be very close to the inner HII region, 
but any conclusion on its properties remains very uncertain. We 
finally point out that it is well known that simple rotational dia- 
grams can overestimate the characteristic rotational temperature, 
sin ce lower transitions typically show higher opacities. 

Mart in-Pintado & Bachillerl (fl992) observed cm-wave inver- 
sion transitions of NH3 up to 7, K-1,1, requiring several hun- 
dred K to be excited. They measured a typical temperature of 
about 150 K in CRL 2688, compatible with our results, but 
without a clear distinction between the nebular components. 
In CRL 6 18, they detect lines mostly in absorption against the 
central free-free continuum source. They were able to identify 
two independent components. One of them, in slow expansion 
and showing a temperature of about 270 K, is very probably 
the slow-expanding inner shell responsible for the narrow fine 
core f ound from Herschel high-7 CO data, see Buiarraba l et al.l 
(l2010h . The other component correspond s to the fast outflow 
also discussed by Buiarrabal et al.l (1201 Oh . and its temperature 
from NH3 data is about 1 5 K. The temperature derived by 
Martin-Pintad o & Bachillerl (1 1992b for the slow component is 
very similar to that of the very inner central region in slow ex- 
pansion studied by means of our CO data; but the NH3 tempera- 
ture of the fast outflow is slightly lo wer than the one found from 
our res ults. We note that the fact that Martin-Pint ado & Bac hiller 
(fl992h detect lines in absorption makes it very difficult to iden- 
tify which part of the fast outflow is actually absorbing the 
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continuum, perhaps the warmest regions of the outflow are not 
placed in front of the very small continuum source and were not 
actually detected in ammonia lines. 

Finally, we mention the observations bv lPardo et alJ d2004l) 
of CRL 6 1 8 in a high number of rotational transitions of HC3N in 
the ground and excited vibrational states. These authors deduce 
the presence of a central, molecule-rich region with 7\ ~ 250- 
300 K. The excitation conditions an d extent of this comp onent 
are very similar to those deduced by Bujarrabal et alJ ( 120101) for 
the inner nebula from the Herschel CO profiles. 

4. Conclusions 

We present observations of FIR/sub-mm molecular lines in ten 
protoplanetary or young planetary nebulae (PPNe, PNe), see 
Table 1, obtained with the heterodyne spectrometer HIFI, on 
board the Herschel Space Telescope. We focused on observ- 
ing CO and H2O rotational lines. We have detected 12 CO and 
13 CO 7=6-5, 7=10-9, and 7=16-15 lines, and we also re- 
port observations of H2O ll,0-l(),l, ll,l-0u,0, 32,1—3 1,2-. 3i,2— 22,1, 

42,2-4i,3, 63,5-62,4, and 734-72,5 lines. 12 CO was detected in 
all the sources. At least one H2O line was detected from all 
but three objects, IRAS 22227+5435, the Red Rectangle, and 
IRAS 17436+5003. In the most intense emitters, we have also 
found other lines within the (wide) receiver bands; the NH3 
7, .£=1,0-0,0 and OH 2 ni /2 3/2-1/2 rotational transitions are 
particularly interesting; we also detected lines of C ls O, H2 ls O, 
vibrationally excited H 2 0, SiO, HCN, and CN. 

The obtained full-band spectra are presented in the Appendix 
(in electronic form only). Selected lines are shown in Figs. 1 to 
14. 

Particularly rich are the spectra of CRL 6 18, CRL 2688, 
NGC7027, and OH 23 1.8+4.2, which have been or will soon 
be discussed in dedicated papers. Seve ral spectra of CRL 61 8 
were presented in a preliminary paper (Bujarra bal et al.ll2010h : 
here we show again those spectra, after careful reanalysis and 
calibration, as well as new recent observations. 

We have performed a simplified interpretation of the CO line 
intensity in terms of general excitation state of the two main 
components of the profiles: the line core and the line wings. It 
is known (Sects. 1, 3) that the line core comes from a slowly 
expanding component that can be identified with the remnant of 
the AGB shell. On the other hand, the line wings come from 
fast bipolar outflows, which are thought to be accelerated and 
heated by shocks during the post-AGB phase. We argue (Sects. 
1 and 3.1) that those FIR/submm lines are the best tool to study 
intermediate-excitation regions, in the temperature ranges 100 
K < 7k < 1000 K, which are not well probed by mm- wave or 
optical/NIR observations. 

In Table 2, we give typical CO intensities of these features, 
as well as intensity ratios between them or between the observed 
lines and the 7=2-1 line (from IRAM 30-m telescope observa- 
tions), both for 12 CO and 13 CO. Significant ratios are represented 
in Figs. 15 to 17, see full description of the parameters in Sect. 3. 
As we can see, in CRL 618 the line wings are particularly intense 
and relatively more intense for the highest-7 transitions, indicat- 
ing particularly high temperatures in t he fas t gas in this source, 
as already found by Bujarra bal et al.l d2010h . In o thers, notably 
in OH 231.8+4.2 and in agreement with results bv lAlcolea et alJ 
d20Tlh . the fast gas seems significantly cooler. 

We also present theoretical predictions of FIR/submm line 
intensities for wide ranges of the kinetic temperature and den- 
sity, performed for a quite general case, see details in Sect. 3.1 
and results in Figs. 18 to 21. From comparison between these 



Table 4. Characteristic temperatures of the fast outflows derived 
in this work, compared with other nebular parameters. 



source 


T k 


flow lifetime 


n lot 


refer. 


CRL 618 


200 K 


< 100 yr 


5 W b cm 3 


2 


CRL 2688 


100 K 


200 yr 


~ 10 6 cnr 3 


1,3 


OH 23 1.8+4.2 


30 K 


1000 yr 


~ 2 10 5 cm" 3 


4 


NGC6302 


40 K 


< 800 yr 


2 10 4 cm" 3 


5 


Frosty Leo 


60 K 


~ 1700 yr 


10 4 cm" 3 


6 


IRAS 17436+5003 


60 K 


~ 500 yr 


~ 10 6 cm" 3 


1,7 



Reference s: 1: this paper (Se ct. 3.11: 2: [ Sanchez Contreras et al] 
J20041); 3: ICox et al] ffijooh: 4: lAlcolea et al] fcOOlh: 5 : from data 
in iDinh-V-Trung et all d2008h: 6: ICastro-Carrizo et al.l d2005h : 7: 
Hoog zaad et al] d2002h . lCasiro-Carrizo et al] d2004h . 



theoretical expectations and our data, we have estimated charac- 
teristic physical conditions of the fast gas in six nebulae, focus- 
ing on the kinetic temperatures. 

We must remember that such a comparison is not straightfor- 
ward, in particular because of the very probable complex struc- 
ture and dynamics of the emitting regions. In fact, in view of 
some general trends in the comparison of our observational and 
theoretical results, we argue (Sect. 3.1) that several components 
may be systematically present and that the densest parts of the 
fast gas may be systematically cooler (see below). 

In the best studied objects, information obtained from pre- 
vious investigations help us to understand the contributions of 
the various nebular components, and our discussion is more re- 
liable. For CRL 618 and OH 23 1.8+4. 2, we can compare our es- 
timates with results from detailed modeling of the FIR/sub-mm 
line emission and we find that both results are quite consistent. 
Therefore, we believe that our estimates of the density and tem- 
perature, although simple, are good approximations to the aver- 
age physical conditions in the fast outflows of the studied nebu- 
lae. We derive a relatively high temperature for the high-velocity 
gas in CRL 618, with typical values 7\ ~ 200 K, probably higher 
in the base of its (relatively compact) bipolar outflow. However 
the fast gas in OH 23 1.8+4. 2 is quite cold, at around 30 K, and 
shows no signs of significantly higher temperatures in any part of 
its very elongated and extended outflow. In CRL 2688, a nebula 
that is also relatively well studied, the characteristic temperature 
of the high-velocity gas is intermediate, 7\ ~ 100 K. Table 4 
shows a summary of the derived characteristic temperatures, as 
well as the time elapsed since the acceleration of bipolar flows 
(or flow lifetime) and an estimate of the characteristic density 
in that component. (In the estimates of the lifetimes and densi- 
ties, we take results from previous studies into account, see ref- 
erences in the Table.) As we can see, there is a clear relation be- 
tween the temperature of the fast gas and the time elapsed since 
its acceleration, with clear extreme values for OH 23 1 .8+4.2 and 
CRL 6 18, and intermediate values for CRL 2688 and other neb- 
ulae. We interpret this result as showing that the gas was ac- 
celerated and heated by the passage of a shock and it is being 
cooled down since then. The density seems also to play a role, 
as expected, since radiative cooling increases with the gas den- 
sity. For instance, the relatively diffuse gas in Frosty Leo and 
NGC 6302 is somewhat warmer than in OH 23 1 .8+4.2, for com- 
parable times. 

We have also discussed briefly other cases of interest: 
NGC 7027 and the Red Rectangle (Sect. 3). NGC 7027 shows 
complex profi les corresponding to it s well known very com- 
plex structure (Nakashima et al. 2010). The molecule-rich shell 
is elongated and very inhomogeneous and shows a composite 
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velocity field. Detailed modeling is absolutely necessary to un- 
derstand the molecular emission from this nebula. The molec- 
ular gas in the Red Rectang le forms an equatorial disk in rota- 
tion tiBuiarrab al et al.ll2005l it is the only PPNe or PNe in which 
gas in rotation has been clearly detected). From a first-approach 
comparison between our profiles in the Red Rectangle and the- 
oretical predictions, we derive a characteristic temperature 7\ > 
150 K (for the gas responsible for the profile maxima), which is 
at least twice the previous estimates from mm-wave data. Due 
to the important role of thermal diffusion in this case (Sect. 3.1), 
the high temperature in the disk must have strong effects on its 
shape and kinematics and can, in particular, explain the presence 
of gas evaporation in the outer parts of the disk. 
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Appendix A: Observational results (in online 
version only) 

In this Appendix we show the full DSB spectra obtained in our 
program. USB and LSB frequencies are indicated in the figures. 
We also show the frequencies of some remarkable molecular 
lines, which does not mean that these lines have been detected in 
the corresponding sources. Depending of the line widths and S/N 
ratio, the spectral resolutions were averaged to channel widths 
ranging between 0.5 and 1.5 km s -1 (except for the narrow lines 
of the Red Rectangle, for which we adopted a resolution of 0.2 
km s ). Selected lines are shown in more detail (and in LSR 
velocity scale) in Figs. 1 to 14. 
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Fig.A.l. Full-band spectra obtained in CRL618. We show r m b vs. upper- and lower-side band frequencies. Frequencies of some 
relevant lines are indicated (which does not mean that they are detected in this source). All frequencies are corrected for Doppler 
shifts due to the relative velocity of the source and therefore correspond to rest frequencies. Continuum has been subtracted in all 
the spectra presented in this paper. 
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Fig. A.2. Full-band spectra obtained in CRL 2688; see caption of Fig. A.l. 
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Fig. A.3. Full-band spectra obtained in NGC 7027; see caption of Fig. A. 1 . 
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Fig. A.4. Full-band spectra obtained in the Red Rectangle; see caption of Fig. A.l. 
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Fig. A.5. Full-band spectra obtained in IRAS 22272+5435; see caption of Fig. A.l. 
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Fig. A.6. Full-band spectra obtained in OH 23 1.8 +4. 2; see caption of Fig. A.l. 
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Fig. A.7. Full-band spectra obtained in the Boomerang Nebula; see caption of Fig. A.l. 
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Fig. A.8. Full-band spectra obtained in Frosty Leo; see caption of Fig. A.l. 
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Fig. A.9. Full-band spectra obtained in IRAS 17436+5003; see caption of Fig. A.l. 
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Fig. A.10. Full-band spectra obtained in NGC 6302; see caption of Fig. A. 1 . 



